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The  principal  investigators,  Greg  Forest  and  Stephen  Bechtel,  have  made  progress 
on  analytical,  numerical,  and  experimental  fronts  in  our  studies  of  non-Newtonian  fluids. 
The  resources  afforded  by  this  grant  have  contributed  significantly  to  these  research 
efforts.  In  addition  to  summer  support  for  the  p.i.'s,  the  AFOSR  has  partially  supported 
Basab  Maulik,  a  postdoc  from  Princeton  University,  Karen  Bolinger,  a  1989  Ph.D. 
graduate  under  the  advisement  of  Forest,  and  Qi  Wang,  a  current  Ph.D.  candidate  of 
Forest's  who  will  graduate  in  1991. 

We  summarize  the  main  results  which  have  been  attained  during  this  granting 
period  and  the  publications  which  have  resulted  that  acknowledge  this  AFOSR  grant. 
The  results  will  be  separated  into  an  analytical  category,  a  numerical  category,  and  an 
experimental  category. 

I.  Analytical  Studies 

We  now  have  produced  two  asymptotic  theories  (pointwise  and  averaged)  for 
slender  free  surface  non-Newtonian  jets.  This  analysis:  (a)  reduces  spatial  dimensions 
from  three  to  one;  (b)  yields  a  seven  parameter  family  of  closed  systems  of  pde's  in  one 
space,  one  time  that  govern  slender  free  non-Newtonian  jets;  (c)  establishes  formal 
consistency  to  all  orders  of  perturbation  for  these  1-D  models;  (d)  produces  explicit 
models  of  a  catastrophic  chan ge-of- type  from  hyperbolic  to  mixed  hyperbolic/elliptic 
type.  These  results  are  contained  in  references  [1-6]. 

These  studies  have  spawned  several  new  ongoing  projects  related  to:  numerical 
loss-of-evolution,  change-of-type  in  the  full  3-D  Maxwell  models,  physical  validity  and 
consequences  of  a  mathematical  change-of-type,  numerical  detection  and  consequences 
of  change-of-type,  physically  meaningful  regularizations  that  produce  well-posedness  in 
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an  ill-posed  regime,  the  analysis  of  convergent  and  stable  numerical  codes,  and  molecular 
dynamics  simulations  of  dilute  polymer  solutions.  These  current  projects  involve 
collaborations  with  M.  Gunzburger,  VPI,  James  Turner,  Hampton  University,  and  C. 
Jayaprakash,  Department  of  Physics,  Ohio  State  University,  as  well  as  two  Ph.D.  students 
of  Forest,  Q.  Wang  and  J.  Cao. 

II.  Numerical  Studies 

We  have  several  interconnected  numerical  studies,  involving  the  p.i.'s,  Maulik, 
Wang,  and  more  recently  Gunzburger  and  Turner.  Wang  has  focussed  on  the  numerical 
simulation  of  specific  asymptotic  model  systems.  These  codes  compute  well-posed 
boundary  value  problems  for  steady  solutions  based  on  the  characteristics  of  the  pde 
systems.  Steady  states  are  computed,  linearized  stability  of  steady  states  and  effects  of 
boundary  conditions  are  computed,  and  the  full  time-dependent  codes  are  developed  and 
run  in  a  neighborhood  of  steady  states  to  explore  basins  of  stability  and  robustness  with 
respect  to  parameters.  The  systems  of  pde's  have  change -of-type  behavior,  and  there  are 
nontrivial  coding  issues  related  to  this.  It  is  here  that  mathematically  deduced  physical 
regularization  of  a  change-of-type  is  critical.  We  know  that  in  our  1-D  model  problems 
the  retardation  effects  regularize  the  elliptic  change  oi- type  from  high  Weissenberg 
numbers,  but  only  in  a  neighborhood  of  constant  solution.,.  The  development  of  a  stable 
code  for  the  regularized  equations  is  in  progress.  This  project  is  prototypical  of  the  3-D 
loss-of-evolution  pathologies. 

Maulik  has  numerically  studied  laminar  transition  in  2-D  and  3-D  Maxwell  fluids. 
These  studies  involve  massive  codes  and  the  results  are  now  beginning  to  get  published. 
We  anticipate  comparison  of  these  macroscopic  simulations  with  molecular  dynamics 
simulations  of  dilute  polymer  solutions  by  C.  Jayaprakash,  Department  of  Physics, 
O.S.U.  The  goal  of  this  project  is  to  model  realistic  dynamics  of  the  large  polymer 
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molecules  and  ultimately  infer  corrections  to  Maxwell  fluid  models  which  are  connected 
to  realistic  polymer  dynamics. 

III.  Experimental  Studies 

Experiments  have  been  performed  involving  free  surface  jets  of  three  very 
different  types  of  fluids:  (1)  water,  (2)  aqueous  biopoljmer  solutions  in  three  different 
concentrations  (Xanthan  gum),  and  (3)  an  aqueous  surfactant  solution 
(Ethoquad/sodium  salicylate).  The  biopolymer  and  surfactant  solutions  are  viscoelastic. 

The  experiments  imply  that: 

1)  robust,  slender  jet  profiles  exist,  so  that  the  assumptions  of  our  asymptotics  are 
supported; 

2)  noticeable  effects  due  to  material  (viscoelastic)  properties  are  documented  in 
otherwise  identical  experiments. 

An  open  problem  in  rheology:  How  does  one  measure  relaxation  time  Aj  in 
elongational  flows?  (Measurement  of  elongational  response  in  fluids  has  proven  difficult 
because  regions  of  elongational  flow  are  in  general  difficult  to  isolate  and  involve  free 
surfaces,  which  complicates  the  analysis.  Elongational  response  cannot  be  inferred  from 
the  more  easily  measured  shear  response  of  a  fluid.)  Our  asymptotic  models  for  free  jet 
flow  (which  is  elongational)  predict  experimentally  measurable  effects  as  functions  of  A  ^ 
for  other  parameters  fixed.  This  poses  an  intriguing  inverse  problem  from  our 
analysis:  we  infer  Aj  by  comparing  experimental  properties  of  the  free  jet  flow  (e.g.,  the 
wavelength  and  amplitude  of  the  observed  nonlinear  spatial  oscillation  in  the  steady 
profile  of  a  jet  issuing  from  an  elliptical  orifice,  and  the  growth  or  decay  of  these 
quantities  in  successive  cycles  down  the  jet)  with  the  corresponding  value  of  Aj  that 
produces  the  same  properties  in  the  1-D  model. 


4 


REFERENCES 


1.  S.E.  Bechtel,  M.G.  Forest,  D.D.  Holm,  and  K.J.  Lin,  "1-D  Closure  Models  for  3-D 
Incompressible  Viscoelastic  Free  Jets:  von  Karman  Flow  Geometry  and  Elliptical 
Cross  Section,"  Journal  of  Fluid  Mechanics,  Vol.  196,  November  1988,  pp.  241-262. 

2.  S.E.  Bechtel,  M.G.  Forest,  and  K.J.  Lin,  "1-D  Closure  Models  for  3-D  Incompressible 
Viscoelastic  Free  Jets:  Higher  Order  Flow  Perturbations,"  Journal  of  Non-Newtonian 
Fluid  Mechanics,  to  appear. 

3.  S.E.  Bechtel,  M.G.  Forest,  K.J.  Lin,  and  D.D.  Holm,  "1-D  Closure  Models  for 
Slender  3-D  Viscoelastic  Free  Jets  with  Elliptical  Cross  Section,"  Proceedings  of  the 
AIAA/ASME/SIAM/APS  1st  National  Fluid  Dynamics  Conference,  July  1988, 
Cincinnati,  Ohio,  Part  3,  pp.  2062-2069. 

4.  S.E.  Bechtel,  K.D.  Bolinger,  and  M.G.  Forest,  "On  the  Validity  of  Strongly 
Viscoelastic  Slender  Jet  Models,"  submitted  to  the  Journal  of  Rheology. 

5.  M.G.  Forest  and  Q.  Wang,  "Change-of-Type  Behavior  in  Viscoelastic  Slender  Jet 
Models,"  submitted  to  Theoretical  and  Computational  Fluid  Dynamics. 

6.  S.E.  Bechtel,  K.D.  Bolinger,  J.  Cao,  and  M.G.  Forest,  "Pointwise  Closure  Models  for 
Slender  Viscoelastic  Free  Jets  with  Torsion,"  to  be  submitted  to  the  Physics  of 
Fluids  A. 

7.  S.E.  Bechtel,  "The  Oscillation  of  Slender  Elliptical  Inviscid  and  Newtonian 
Jets:  Effects  of  Surface  Tension,  Inertia,  Viscosity,  and  Gravity,"  AS  ME  Journal  of 
Applied  Mechanics,  Vol.  56,  December  1989,  pp.  968-974. 

8.  B.K.  Maulik,  "Linear  Instability  of  Laminar  Flow  for  the  Oldroyd-B  Fluid  in  a  Planar 
Channel,"  Journal  of  Non-Newtonian  Fluid  Mechanics,  to  appear. 


5 


PRESENTATIONS 

"On  the  Behavior  of  Viscoelastic  Free  Jets  with  Elliptical  Cross  Section,"  Applied 
Mechanics  and  Engineering  Sciences  Conference,  University  of  California,  Berkeley, 
June  1988;  abstract  in  Section  7  of  the  proceedings. 

"1-D  Closure  Models  for  Slender  3-D  Viscoelastic  Free  Jets  with  Elliptical  Cross 
Section,"  AIAA/ASME/SIAM/APS  1st  National  Fluid  Dynamics  Congress,  Cincinnati, 
Ohio,  July  1988;  paper  in  the  proceedings  on  pp.  2062-2069. 

"1-D  Closure  Models  for  Newtonian  and  Viscoelastic  Free  Jets:  Lowest  Order  Equations 
and  Higher  Order  Corrections,"  41st  Annual  Meeting  of  the  American  Physical  Society, 
Division  of  Fluid  Dynamics,  Buffalo,  New  York,  November  1988;  abstract  on  page  2245 
of  the  program. 

"The  Effects  of  Viscoelastic  Properties  on  the  Oscillation  of  Elliptical  Free  Jets,"  60th 
Annual  Meeting  of  the  Society  of  Rheology,  Gainesville,  Florida,  February  1989; 
abstract  on  page  40  of  the  program. 

"Asymptotic  Closure  Models  for  Non-Newtonian  Fluid  Jets,"  invited  talk  at  the  Applied 
Mathematics  Colloquium,  University  of  Colorado,  Boulder,  September  1989. 

"Experimental  and  Theoretical  Studies  of  the  Behavior  of  Viscoelastic  Free  Jets  from 
Elliptical  Orifices,"  61st  Annual  Meeting  of  the  Society  of  Rheology,  Montreal,  Canada, 
October  1989. 

"The  Oscillation  of  Elliptical  Inviscid  and  Newtonian  Jets:  Effects  of  Surface  Tension, 
Inertia,  Viscosity,  and  Gravity,"  Winter  Annual  Meeting  of  the  ASME,  San  Francisco, 
California,  December  1989. 

"Nonlinear  Oscillations  in  Newtonian  and  Viscoelastic  Free  Surface  Jets,"  invited  talks  at 
the  Ohio  State  University  Engineering  Mechanics  Seminar,  Columbus,  February  1990, 
and  the  University  of  Michigan  Mechanical  Sciences  Division  Seminar,  Ann  Arbor, 
March  1990. 

"Observations  and  Analysis  of  Newtonian  and  Non-Newtonian  Free  Surface  Jets," 
invited  talk  at  the  University  of  Rochester  Mechanical  Engineering  Colloquium, 
Rochester,  New  York,  March  1990. 

"Newtonian  and  Viscoelastic  Free  Jet  Behavior,"  to  be  presented  at  the  Eleventh  U.S. 
National  Congress  of  Applied  Mechanics,  Tucson,  Arizona,  May  21-25,  1990. 

"On  the  Validity  of  Viscoelastic  Slender  Jet  Models,"  to  be  presented  at  the  62nd  Annual 
Meeting  of  the  Society  of  Rheology,  Santa  Fe,  New  Mexico,  October  21-25,  1990. 


